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I. INTRODUCTION 

Multiple scattering occurs when waves propagate in media with scatterers, leading to 
many interesting phenomena such as the bandgaps in periodically structured media and 
the Anderson localization in disordered media[l|. Within a bandgap, waves are evanescent; 
when localized, they remain confined spatially near the transmission point until dissipated. 
The phenomenon of bandgaps and localization has been both extensively and intensively 
studied for electronic, electromagnetic, and acoustic systems, as summarized in Ref. j^. 

Propagation and scattering of gravity waves over topographical bottoms has also been 
and continues to be a subject of much research. From practical perspectives, it is essential 
to many important ocean engineering problems such as designing underwater structures in 
minimizing the impact of water waves on banks. From the point of fundamental research, the 
system of water waves may also play an important role. Since water waves are a macroscopic 
phenomenon, they may be monitored and recorded in a laboratory scale. In this way, many 
significant phenomena of microscopic scales may be demonstrated with water waves. This 
would be particularly useful in facilitating the understanding of abstract concepts which 
may have been initiated for quantum waves. A great amount of papers and monographs 
has been published on water waves over various topographical bottoms 13, 14 S la ID, la Isl 

um\~~ 

textbooks 



15 



16 



13, 14 1 . A comprehensive reference on the topic can be found in three excellent 



id. 



The concept of bandgaps and Anderson localization has been extended to the the study 
of the propagation of water surface waves over topographical bottoms. In 1983, Guazzelli 
et al. |18| suggested that the phenomenon of Anderson localization could be observed 
on shallow water waves, when the bottom has random structures. Later, Devillard et al. 
reconsidered water wave localization inside a channel with a random bottom in a framework 
of the potential theory 19] . They computed the localization length for various cases. The 
experimental observation of water wave localization has been subsequently suggested by 
Belzone et al. \2(^ . 

When the topographical bottoms are periodically structured, the propagation of water 
surface waves will be modulated accordingly. According to Bloch theorem, waves in a 
periodic medium, now termed as Bloch waves, can be expressed in terms of the product of a 
plane wave and a periodic function which is has the periodicity of the medium. Therefore, the 



waves will exhibit the properties of both plane wave propagation and periodic modulation. 



Indeed, a recent experiment 2l| has used gravity waves to illustrate the phenomenon of 
Bloch wave over a two dimensional periodic bottom. This pioneering experiment has made 
it possible that the abstract concept be presented in an unprecedentedly clear manner. The 



22|. 



experimental results have also been matched by a theoretical analysis in Ref. 

The experimental advanturesJ20|, |2l| pave a new avenue for investigating the phenomena 
of Anderson localization in disordered media and wave bandgaps in periodic structures. In 
particular, it has been recognized that the two dimension is a marginal dimension, many 
peculiar wave phenomena could be expected J2j]. Since gravity waves are naturally a two 
dimensional phenomenon, it may be expected that water surface waves can play a role in 
demonstrating some of the peculiar wave phenomena. 

Motivated by the developments, we wish to further consider the propagation and localiza- 
tion properties of water surface waves. At the first stage, we will consider water waves over 
one dimensional uneven bottoms. We present a theoretical analysis of the previous experi- 
mental results |2nl|. The formulation in Ref. [2^ will be used for the purpose. The comparison 
between the experimental and theoretical results, in return, provides a verification of the 
theory. We will study the band structure of periodic cases, the effect of randomness on wave 
propagation, the relation between the bandgaps and the localization, and the amplitude or 
energy distribution over the structured bottoms. The dependence on parameters, such as 
the frequency, water depth, the variations of the height and width of the obstacle steps, will 
be examined in detail. Although the experiment to be analyzed here was done nearly twenty 
years ago, to the best of our knowledge, however, there are no further experiments which 
have been done on water waves in the context of the localization effects. Even the existing 
limited experimental results have not been thoroughly analyzed. The present paper is to 
bridge the gaps, with the hope that further experimental investigation may be arranged. 

The paper will be constructed as follows. In the next section, we will present the formu- 
lation and parametrization for the problem. The results and discussion will be presented in 
section III, followed by a summary in section IV. 



II. THE GENERAL FORMULATION 



A theory of water wave propa. 



posed and developed in Refs. 
some experimental data 



m 



^ion over step-mounted bottoms has been recently pro- 

23|. This formulation has been used earlier in interpreting 

22j |. While the details can be referred to in Ref. J23[, here we only 

present the final equation. After the Fourier transformation, the equation describing the 

wave of frequency uj over topographical bottoms is 

^ (^^^) + ^ = 0' (1) 

where k satisfies 

cj^ = gk{r)ta.iah{k{r)h{r)), (2) 

where t] is the surface displacement, g is the gravity acceleration constant, and h is the 
depth from the surface. For a fixed frequency, the variation of the wavenumber k with the 
topographical bottom is determined by the depth function h. 

From Eq. (^, we have the conditions linking domains with different depths as follows: 
both r] and ^^^^ — -i] = -j^t] are continuous across the boundary. 

A. Application to one-dimensional situations 

1. A single step 

First, consider a step with width d, and a wave is propagating along the x direction. The 
conceptual layout is in Fig. ^a). We use the standard transfer matrix method to solve for 
the wave transmission across the step. 

The waves on the left, within, and on the right side of the step can be generally rewritten 
as 

riR = Arc''^^^ + BRe-'^^\ (3) 

The subscripts L, M, R represent the quantities on the left side, in the middle, and on the 
right side of the step respectively. 

The boundary conditions lead to the following equations: 



and 
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(7) 



Km ^ Kr 

In these equations, xl^r stands for the locations of the left and right sides of the step 
respectively, and xr — xl = d. 

The first set of boundary equations gives the matrix relation 
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From Eqs. (jH)) and ()11|) . we obtain the following solution in the transfer matrix form 

Ar 




T, 



LR 



B 



(12) 



R 



with 

Tlr = TlmTmr- (13) 

Eq. (fT^ relates the waves on the left to the right side of the step. 

2. The case of N steps 

Now we consider A^ steps in a unform medium of wave number k. The illustration is in 
Fig. n^b). The step widths are di and the water depths over the steps are h^. The wave 
number over the step is denoted by /cj {i = 1, . . . , N). They satisfy the following relations 
respectively, 

oj"^ = ghtanh-ikh), oj"^ = ghitaiih.{kihi). (14) 

Clearly, the coefficients on the most left hand region is related to the most right hand 
region by 

h\^nN)(^A (15) 

with 

N 

T(iV) = nT, (16) 

i=l 

The matrix Tj is the transfer matrix for the z-th step and will be given below. 

Let us consider a unit plane wave propagation along the x direction, and explore the 
reflection and transmission properties. In this case, clearly we have 

Al = 1, Br = 0. (17) 

Br = is the common radiation condition. Thus from Eq. (fT3j) we arrive at the solutions 

The subscripts denote the corresponding matrix elements. 
The transmission and reflection coefficients are defined as 

T = \Ar{N)\\ R=1-T. (19) 

Now we construct the T matrix for each step. In the current case, we have 

k k- 

9LM{i) = -r, 9MR{i) = -77, (20) 

/Cy fv 



and 

kh = k, ku = ki, kn = k. (21) 



We denote Qsi = t'- Therefore 
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B. Simulation setup 

1. Non-dimensional parametrization 

Consider an infinite periodic array of the steps, as shown in Fig. ^ The lattice constant 
is L. For random arrays, L refers to the average distance between two adjacent steps. 
The dispersion relation is 

uj^ = gktajih.{kh). 

This can be rewritten as 

— = (A;L)tanh {kL)-\ , 
uJq \ LI 



with 
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Therefore in all later computations, the length can be scaled by L, the frequency by Uq, the 
wavenumber by kL. 

The wavenumbers in the medium and within the steps are given by (at the same fre- 
quency) 

^2 f h 

-^ = (A;L)tanh (kL)- 



2 - (A;L)tanh( [kL)-] , (23) 
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and the transfer matrix of the 


i-th step is 
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(25) 



where Xj^^, is the coordinate of the left side of the z-th step. 

2. Band structure for periodic cases 

For the periodically arranged steps with di = d, and hi = hi, the band structure can be 
solved. According to Bloch theorem, the water surface displacement r] can be expressed as 

V{x) = ^(a;)e^^^ (26) 

where K is the Bloch wavenumber, and ^(x) is a periodic function modulated by the peri- 
odicity of the structure, i. e. ^(x + L) = C,{x). The relation between K and the frequency lo 
can be obtained by taking Eq. (j^ into Eq. (fT^. 

We can derive an equation determining the band structure in the periodic case, 

cos{KL) = cos{kiL{d/L) cos{kL{L - d)/L) - cosh(20 sm{kiL{d/L)) sm{kL{L - d)/L), 

(27) 
where 

^ = ln(g), with u'^ = gkita.nh.{khi), q^ 



,2 L +„„l,/'L^ N „2 _ ^1 

" ~ k' 



3. Random situations 

There are a number of ways to introduce the randomness. (1) Variation in the height 
of the steps: With the fixed widths and positions of the steps, the height of the steps can 
be varied in a controlled way. For example, the height of the steps can be varied randomly 
between [Hq — AH, Hq + AH]. (2) Positional disorders: Initially, the steps can be arranged 
in a lattice form. Then allow each step to move randomly around its initial position. The 
allowing range for movement can be controlled and denotes the level of randomness. The 

8 



extreme case is completely randomness. (3) Width randomness: We can also introduce the 
randomness for the widths of the steps. In the simulation, we will consider the randomness 
introduced in the experiment [20| . 

When randomness is introduced, a few definitions are in order. The most important 
quantity is the Lyapounov exponent 7. Its definition is 



where 



7=lim(7jv), (28) 



7^ = -^ln(|A^(iV)|). 



Here |Ak(A^)P is the transmission coefficient for a system with A^ random steps, referring 
to Eq. P|). 

\MN)\'- ^ 



|Tn(iV)r 

The symbol (■) denotes the average over the random configuration. The inverse of the 
Lyapounov exponent characterizes the localization length, i. e. ^ = 7^^. 

III. THE RESULTS AND DISCUSSION 

The systems are from the previous experiment |2Cll|. That is, the bottoms are mounted with 
a series of steps and these steps are either regularly or randomly but on-average regularly 
placed on the bottoms. Three cases are considered and are illustrated by Fig. |21 In the 
Bed P case, the averaged water depth is Hq, the periodicity is 2Lo, and the step variation 
is fixed at aH. In the Bed RS case, the averaged water depth is Hq, the step variation is 
fixed at aH, the separation between steps is uniformly distributed with [cIq — Ad, do + Ad]. 
In the Bed R case, both the height and the separation between the steps are allowed to vary 
randomly, but within the ranges [Hq — AH, Hq + AH] and [Lq — AL, Lq + AL] respectively. 

The experimental setups have been described in Section 2 of Ref. |20j. We briefiy repeat 
here. The experiments were carried out in a glass-walled wave tank with length 4m and 
width 0.39m. A bottom composed of periodic or random steps was built into a fiat bottom 
with the mean water depth Hq. The different bottoms varied only along the tank so that, 
apart from weak edge effects, the propagation of waves is considered to be one-dimensional. 
The resolution of the water depth is estimated at about 0.2mm. 



A. Band structure and transmission 

First, we consider the first case in the experiment: the periodic case, i. e. the Bed P case. 
For this case, the band structure and the transmission are computed for two water depths. 
In both cases, the width of the steps is Lq = 4.1cm; therefore the periodicity is 8.2cm. The 
results are shown in Fig. 01 From the band structure results in (al) and (bl), we observe that 
for the small water depth {Hq = 1.75cm), there are two bandgaps in the frequency range 
measured in the experiment, while in the deeper case {Hq = 3cm) there is one bandgap. The 
locations of the band gaps match the inhibited transmission regimes. The width of the gap 
and the inhibition effect tend to decrease with frequency as shown by (al) and (a2). This is 
understandable. In the high frequency limit, i. e. when kh » 1, the dispersion relation in 
Eq. (J2I) approaches u"^ = gk. Therefore the importance of the bottom structure will decrease 
with increasing frequency. The results in Fig. El will help us comprehend the later results. 

B. Reflection coefficient 



In the experiment [20|] , the reflection coefficients are measured for the three cases: Bed 
P, Bed R and Bed RS cases. Two average water depths are considered: i^o = 1-75 and 3 
cm. We have considered all the cases, and applied the formulation in Section II to obtain 
corresponding results. In Fig. 01 we present our theoretical results. For the convenience of 
the reader and as a comparison, we also re-plot the experimental results on the same figure 
(left panel). We have take into account two random configuration numbers in the simulation: 
one is 5 random configuration, i. e. the middle panel, which is take as the same as in the 
experiment; the other in the right panel is more than ten thousand random configurations to 
ensure the stability of the averaging results. All the parameters are repeated from Ref . J2i| . 

Fig. m^al), (a2) and (a3) compare the results for the Bed RS and Bed P cases with 
averaged water depth Ho = 1.75 and step width Lq = 4.1cm. For both cases, the ratio 
aH /Hq is fixed at 0.43, i. e. the there is no variation in the step heights. In the Bed RS 
case, the disorder is introduced to the separation between steps, that is, the separation is 
randomly chosen with uniform distribution within [2cm, 8cm] or [(io— 3, (io+3] with d^ = 5cm. 
In the simulation, the total number of steps is 58. We have taken two numbers of random 
configuration in the simulation. One is five (a2), which complies with the experiment, 
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and the other (a3) is ten thousand time, so to ensure the stabihty of the averaging. The 
experimental data are shown in (al). The comparison of (al), (a2) and (a3) indicates the 
following. Overall speaking, the theoretical results capture well the qualitative features 
observed experimentally, and agree to certain extend with the experimental results. 

First we consider the Bed P case. (1) The theoretically predicted positions of the reflection 
peaks agree well with the experimental observation in the Bed P case. These positions also 
coincide the bandgaps from the band structure computation in Fig. ^ (2) In the Bed P case, 
the reflection coefficient reaches its maximum value of one within the bandgaps as expected, 
while the experimental values are always smaller than one for the frequency range considered. 
A possible reason for this discrepancy may be that in the theoretical simulation, we did 
not take into account possible dissipation effects caused by such as viscosity and thermal 
exchange; some effects have been discussed in Ref. |2fl]. (3) The theoretical width of the first 
refiection peak in the Bed P case matches well that observed, but the theoretical width of the 
second refiection peak at about 4 Hz is narrower than that from the experiment. In fact, the 
experimentally measured widths of the two refiection peaks are more or less the same. Since 
the effects of the periodic bottom diminish with increasing frequency as discussed above, 
we may conclude that there are other effects which could broaden the reflection peak at 4 
Hz, and these effects may include those from the dissipation, non-linearity, and evanescent 
mode leakage. These effects have not been considered in the current theoryJ23|. 

Now we consider the Bed RS case. (1) Again, except for the peak values in the reflection, 
the theoretically results reproduce the experimental observation reasonably well in general, 
particularly at the strong reflection located at about 2Hz. (2) Different from the Bed P case, 
the width of the flrst reflection peak at 2Hz is wider in theory than in the experiment. (3) 
In both the theory and experiment, a second reflection peak is noticed within 4 and 5 Hz. 
(4) An obvious difference between the theory and the experiment is at the low frequency 
around 1.2Hz: A strong sharp reflection peak appears at 1.2Hz in the experiment, but absent 
in the theory. This experimental observation differs from previous observation in acoustic 
or optical systems [2] • In acoustic or optical systems, the disorder effect decreases as the 
frequency decreases. Therefore, waves tend to diffuse away, leading to weaker reflections at 
low frequencies. In fact, the result of the reflection measurement is also in disagreement 
with the localization measurement shown in Fig. 15 of Ref. [20| where it is shown that the 
localization length at 1.2Hz is even longer than at 1.5Hz at which the reflection is small; 
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the longer the locahzation length, the weaker is the reflection. (5) Increasing the number of 
random configuration tend to smooth the curves. 

The comparison between the theoretical and experimental reflection results for the Bed 
P and Bed R cases with Hq = 3cm is shown by Fig. 0I^bl), (b2) and (b3). The parameters 
are as follows. (1) Bed P: aH/Ho = 0.25, Lq = 4.1cm; (2) Bed R: the separation between 
steps varies completely randomly within Lq ± AL with AL = 2cm, and the height of the 
steps varies uniformly within Hq ± AH with AH = 1.26cm. The number of steps is 58. In 
(b2), five random configurations are used for averaging, and in (b3) ten thousand random 
configurations are used to ensure the stability of the averaging. In the Bed P case, except at 
the reflection peak, the theoretical results reproduce very well the experimental observation. 
In the Bed R case, the theoretical results also match that from the experiment in both the 
qualitative structure and the magnitude, referring to (bl) and (b2). The existing deviation 
may result from the insufficient random average. 

The comparison between the theoretical and experimental reflection results for the Bed 
P and Bed R cases with Hq = 1.75cm is shown by Fig.EJ^cl), (c2) and (c3). The parameters 
are as follows. (1) Bed P: cxH/Hq = 0.43, Lq = 4.1cm; (2) Bed R: the separation between 
steps varies completely randomly within Lq ± AL with AL = 2cm, and the height of the 
steps varies uniformly within Hq ± AH with AH = 1.24cm. The number of steps is 58. In 
(c2), flve random conflgurations are used for averaging, and in (c3) ten thousand random 
configurations are used to ensure the stability of the averaging. The Bed P case has been 
discussed in the above. In the Bed R case, the general features of the experimental and 
theoretical results seem to be agreeable with each other. The predicted reflection curve 
starts to match qualitatively the experimental data from about 3 Hz. The discrepancy at 
low frequencies is, again, noticeable. 

C. Localization length 



In the experiment [20| , the localization length is extracted from the measurement of the 
total wave amplitude attenuations. In the simulation, the localization length is obtained 
from the inverse of the Lyapounov exponent given in Eq. ()28|1 . Here the Bed R case is 
considered and the parameters are: Hq = 1.75cm, Lq = 4.1cm, and the height of the steps 
and the separation between steps vary randomly within the ranges [Hq — AH, Hq + AH] 
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and [Lq — AL, Lq + AL] respectively; here AH = 1.2425cni, and AL = 2cni. Ten thousand 
steps and ten thousand random configurations have been used in the simulation to ensure 
the stability of the numerical results. 

The numerical and experimental results are shown in Fig. El Here the localization length 
is plotted against the frequency. The results from Ref. [20| are shown in the inserted box. A 
few observations are in order. (1) In Ref. '2(|, the authors have used a potential formulation 
to obtain the localization length, denoted by the solid length in the inserted box. By eye- 
inspection, we see that the present numerical results agree remarkably well with the results 
from the potential theory, thus providing another support for the present relatively simple 
theory, stemmed from Ref. l23l|. (2) The numerical results also agree with the averaged 
experimental data in the vicinity of the frequency 2Hz. (3) There is a huge fiuctuation in 
the experiment results. From our simulation, we think that such a significant deviation 
is due to the insufficient average numbers, an obvious limitation on any experiment. This 
is particularly an important factor when the localization length is long. Nevertheless, the 
agreements shown in Fig. El is encouraging. 

D. Behavior of the wave amphtude along the random bed 

In the experiment, the variation of the wave amplitude along the random bed is also 
measured. Both Bed RS and Bed R cases are considered. The parameters used in the 
experiment [20| are summarized as follows. In the Bed RS case, Hq = 1.75cm, aH/H = 0.43, 
Lq = 4.1cm, and the separation between steps varies randomly in the range of [2cm, 8cm]. 
In the Bed R case, Hq = 1.75cm, Lq = 4.1cm, and the height of the steps and the separation 
between steps vary randomly within the ranges [Hq — AH, Hq + AH] and [Lq — AL, Lq + AL] 
respectively; here AH = 1.2425cm, and AL = 2cm. Four different frequencies have been 
measured and simulated. 

The experimental and simulation results for a given random realization of the random 
beds are presented in Fig. JHl It is shown that the theoretical results match remarkably well 
the experimental results. It is shown that the waves do not decay monotonically along the 
random bottom (without averaging), due to the manifestation of resonant modes of the beds. 
The resonances are sensitive to the frequency variation. We also found that the occurrence 
of the resonances are sensitive to the random configuration. 
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We have further computed the averaged variation of the wave amphtude along the random 
bed for a sufficiently large number of random configurations. We found that though smeared 
out a little by the averaging, the resonance feature remains for spatial points near the 
transmission and tends to diminish for large travelling paths. And the averaged amplitude 
decays exponentially with increasing travelling distances. As an example, in Fig. [7| we 
illustrate these by the results of the Bed RS case with f = 1.6Hz. The results in Fig. [7| 
also indicate that the exponential decay rate, associated with the localization length, may 
not be accurately obtained from measurements done on insufficiently long samples, as the 
fluctuation can be quite signiflcant for small sample sizes. 

IV. SUMMARY 

In summary, we have considered the propagation of water surface waves over topograph- 
ical bottoms. A transfer method has been developed to compute the wave fleld along the 
propagating path, the transmission and reflection coefficients. The localization effects due 
to disordered bottom structures are also considered. The theory has been applied to ana- 
lyze the existing experimental results. Some agreements and discrepancies are discovered 
and discussed. It is pointed out that more detailed experiments may be helpful in not only 
identifying the peculiar localization phenomenon, but in helping improve theories for water 
wave propagation over rough bottoms. 
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FIG. 1: Conceptual layouts. 
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FIG. 2: Situations considered in the paper, adopted from Fig. 2 of Ref. "20"]. (a) Bed P case: in 
this case, the steps are mounted periodically with lattice constant 2Lq; the variation of the steps 
aH is fixed; (b) Bed RS case: in the case, the steps are allowed to move randomly from their 
initial periodic positions, as set in the Bed P case - the allowed range is denoted as itAL and the 
variation of the steps aH is fixed; (c) Bed R: in this case, both the heights and the widths are 
allowed to vary randomly from their initial values in the Bed P case within the ranges [Hq — AH, 
Ho + AH], and [Lq - ALq, Lq + ALq]. 
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FIG. 3: Band structure and transmission for the Bed P case in two situations, referring to Fig.[21^a): 
(al) and (a2) the average water depth is Hq = 1.75cm, the height variation is aH/Ho = 0.43; (bl) 
and (b2) the average water depth is Hq = 3cm, the height variation is ctH/Hq = 0.25. The left 
and right panel show the band structure and transmission results respectively. The transmission 
is presented in the log scale for 100 steps. 
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FIG. 4: Reflection versus frequency for the Bed RS and P cases with three average water depths. 
Left panel: The results from the experiment |2g]. Middle panel: the theoretical results with the 
average over five random configurations. Right panel: the theoretical results with the average over 
ten thousand random configurations, so to make sure the stability of the average. 
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FIG. 5: Localization length versus frequency for the Bed R case. The experimental results are 
shown in the inset figure. The five black dots denote the results from an averaging over five 
random configurations. 
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FIG. 6: Variation of the amplitude of wave elevation along the wave tank for the Bed RS and Bed 
R cases for different frequencies. The experimental J20| and numerical results are shown on the 
left and right panels respectively. 
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FIG. 7: The averaged variation of the amphtude of wave elevation along the wave tank for the Bed 
RS from Fig. El with / = 1.6Hz. To show the behavior near the transmission site, the results are 
plotted in two length scales: (al) up to 2m; and (a2 

up to 20m. 
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